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We evaluate in this work the hydrodynamic transport coefficients of a granular bi- 
nary mixture in d dimensions. In order to eliminate the observed disagreement (for 
strong dissipation) between computer simulations and previously calculated theoreti- 
cal transport coefficients for a monocomponent gas, we obtain explicit expressions of the 
seven Navier-Stokes transport coefficients with the use of a new Sonine approach in the 
Chapman-Enskog theory. Our new approach consists in replacing, where appropriate in 
the Chapman-Enskog procedure, the Maxwcll-Boltzmann distribution weight function 
(used in the standard first Sonine approximation) by the homogeneous cooling state dis- 
tribution for each species. The rationale for doing this lies in the fact that, as it is well 
known, the non-Maxwellian contributions to the distribution function of the granular 
mixture become more important in the range of strong dissipation we are interested in. 
The form of the transport coefficients is quite common in both standard and modified 
Sonine approximations, the distinction appearing in the explicit form of the different col- 
lision frequencies associated with the transport coefficients. Additionally, we numerically 
solve by means of the direct simulation Monte Carlo method the inelastic Boltzmann 
equation to get the diffusion and the shear viscosity coefficients for two and three di- 
mensions. As in the case of a monocomponent gas, the modified Sonine approximation 
improves the estimates of the standard one, showing again the reliability of this method 
at strong values of dissipation. 



1. Introduction 

The success of the kinetic theory tools in describing granular gases has been widely rec- 
ognized (Goldhirsch 2003 Brilliantov & Poschcl 2004). In particular, the Navier-Stokes 
(NS) constitutive equations for the stress tensor and the heat flux of a monocomponent 
gas have been derived ( Sela fc Goldhirsch|1998 Brey et aTp 998) by solving the inelastic 
Boltzmann equation by means of the Chapman-Enskog (CE) method to first order in 
the spatial gradients (Chapman & Cowling 1970). As in the elastic case, the NS trans- 



port coefficients are given in terms of the solutions of linear integral equations (Brey 
et «L|1 998 ) , which can be approximately solved to get the analytical dependence of these 
coefficients on dissipation. The standard method consists in approaching the solutions 
to these integral equations by the Maxwell-Boltzmann distribution /m times truncated 
Sonine polynomial expansions. For simplicity, usually only the lowest Sonine polynomial 
(first Sonine approximation) is retained (we will refer henceforth to this procedure as the 
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standard Soninc approximation). In spite of this approximation, the results obtained from 
this approach (which formally applies for all values of the coefficient of restitution, see 
Brey et al]|1998[) compare very well with Monte Carlo simulation s (Brey, Ruiz-Montero 



fc Cubero|i999t|Brey, Ruiz-Montero, Cubero fe Garcfa-Rojo|2000| [Lutsko, Brey fe Dufty 



2002 



Garzo fe Montanero||2002 1 for mild degrees of inelasticity, namely, coefficients of 
restitution a that are larger than about 0.7. Specifically, while the self-diffusion and 
shear viscosity coefficients are well estimated by the first Sonine approximation (even 
for values of a smaller than 0.7), the transport coefficients associated with the heat flux 
present significant discrepancies with computer simulations ( Brey fe Ruiz-Montero|200~4 



Brey et at. 2005a[ ) for high dissipation (a < 0.7). This fact has motivated the search 
for alternative methods which provide accurate estimates for the NS transport coeffi- 
cients in the complete range of values of a. Thus, Noskowicz et al. (2007) have used 
Borel resummation to obtain the distribution function of the homogeneous cooling state 
(HCS) and the NS transport coefficients. The method involves the solution of a system 
of algebraic equations that must be numerically solved. A different alternative consists 
in a slight modification of the Sonine polynomial approach, assuming that the isotropic 
part of the first order distribution function is mainly g overned by th e HCS distribution 
/(°) rather than by the Maxwellian distribution /m (Lutsko 20051 (and this is what 
we call modified Sonine approximation). Except for this, this modified Sonine approx- 
imation keeps the usual structure of the standard Sonine approximation. Comparison 



with computer simulations ( 


Brey & Ruiz-Montero 


2004 


Brey, Ruiz-Montero, Maynar 


& Garcia de Soria 


2005 a Montanero, Santos & Garzo 


20071 shows that this modified 



approximation significantly improves the accuracy of the NS transport coefficients of the 
single granular gas at strong dissipation (Garzo, Santos & Montanero 2007c), especially 



in the case of the heat flux transport coefficients. The improved agreement in the range 
of high inelasticities is of physical relevance, since it helps to put into proper context the 
discussion vividly maintained in the field (Aranson & Tsimring 2006 Goldhirsch]|2003 1 



about the existence (or not) of scale separation in granular gases. Moreover, the range of 



high inelasticities has growing interest in experimental works ( 


Clerc, Cordero, Dunstan, 


Huff, Mujica, Risso & Varas 2008 Vega Reyes & Urbach 


2008a 1. Thus, it is important 



from a fundamental point of view to understand the origin of the discrepancies between 
theory and simulation for a < 0.7. 

Motivated by the good agreement found between theory and simulation for a mono- 
component granular gas, we extend in this work the modified Sonine approach to the 
case of a granular binary mixture. Early attempts to obtain the NS transport coefficients 
were carried out by means of the CE expansion around Maxwellians at the same temper- 



ature for each species (Jenkins & Mancini 1989 Zamankhan 1995 Arnarson & Willits 



1998 


Willits & Arnarson 


1999 


Serero, Goldhirsch, Noskowicz & Tan 


200G 


). However, 


as confirmed by computer simulations (Montanero & Garz6|2002 Barrat & Trizac|2002 



bot1[2003 Wang, Jin fe~Ma 2003;|Brey, Ruiz-Montero fe Moreno 2005& 2006; Schrotcr. 



Feitosa 



Ulrich, Kreft, Swift fe Swinney|2006[ ) and experiments QWildman fe Parker|2002 
& Menon 20021, the equipartition assumption is not in general valid in granular mix- 
tures since it is only close to being fulfilled in the quasiclastic limit. A kinetic theory for 
granular mixtures at low density which accounts for nonequipartition effects has been 



developed in the past few years by Garzo & Dufty (20021. As in the single gas case, 
the above theory ( Garzo fc Dufty|[2002 l shows that the NS transport coefficients of the 
mixture are given in terms of the solutions of a set of linear integral equations, which are 
approximately solved (Garzo & Dufty 2002 Garzo, Montanero fc Dufty]|2006 Garzo & 
Montanero|2007[ ) by means of the standard Sonine approximation. Furthermore, this the- 
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ory has successfully predicted, as shown by computer simulations (Montanero & Garzo 
2002]|Dahl, Hrenya, Garzo fe Dufty|2002| [Brey, Ruiz-Montero fc Moreno|2005&||Schroter, 



Ulrich, Kreft, Swift & Swinney 2006), a constant ratio between the granular tempera- 



tures of both species. Here, we revisit the theory of Garzo & Dufty (20021 and solve the 
corresponding linear integral equations defining the NS transport coefficients by taking 
the HCS distribution f^ (i = 1, 2) of each species as the weight function. 

As expected, the problem for a binary mixture is much more involved than in the 
monocomponent case since not only the number of transport coefficients is larger but 
these coefficients are also functions of more parameters (masses, sizes, composition and 
three coefficients of restitution) . In spite of these technical difficulties, explicit expressions 
for the seven relevant transport coefficients of the binary mixture (the mutual diffusion 
D, the pressure diffusion D pi the thermal diffusion D' , the shear viscosity 77, the Dufour 
coefficient D", the thermal conductivity A, and the pressure energy coefficient L) have 
been obtained in terms of the parameter space of the problem. The results show that 
the modified Sonine approximation at large inelasticities introduces significant, moderate 
and slight corrections to the transport coefficients associated with the heat flux (D", A, 
L), stress tensor (77) and mass flux (D, D p , D'), respectively. In order to assess the 
degree of accuracy of the modified Sonine approximation, we have also performed Monte 
Carlo simulations for the mutual diffusion and the shear viscosity coefficients and use 



available simulation data for the heat flux coefficients in the monodisperse gas case ( Brey 
fc Ruiz-Montero|2004 Brey, Ruiz-Montero, Maynar fc Garcia de Soria|2005a| . We show 



that, in the range of strong inelasticity, the modified Sonine approximation provides 
better agreement with simulations in all cases, being this accuracy gain not negligible for 
viscosity, thermal conductivity, pressure energy coefficient and Duffour coefficient. This 
clearly justifies the use of this new approximation in order to obtain accurate expressions 
for the NS transport coefficients of the mixture. 

An important issue is the applicability of the NS equations for granular binary mixtures 
derived here. The expressions of the seven NS transport coefficients are not restricted 
to weak inelasticity and hold for arbitrary values of the degree of dissipation. In fact, 
our modified Sonine approximation turns to be more accurate than the previous one 



( Garzo fc Dufty|2002[) for very dissipative interactions. However, as already pointed out 



by Garzo et al. (2006), the NS hydrodynamic equations themselves may or may not be 
limited with respect to inelasticity, depending on the particular granular flows of real 
materials considered. It is important to recall that the CE method assumes that the 
relative changes of the hydrodynamic fields (partial densities, mean flow velocity and 
granular temperature) over distances of the order of the characteristic mean free path of 
the system are small. In the case of molecular or ordinary gas mixtures, the strength of 
the spatial gradients is controlled by external sources (the initial or boundary conditions, 
driving forces), and for a great variety of experimental conditions of practical interest 
the NS hydrodynamic equations describe adequately the hydrodynamics of the system. 
Unfortunately, the situation in the case of granular gases is more complex, since in most 
problems of practical interest (such as steady states) the boundary conditions cannot 
completely control by themselves the spatial gradients of the hydrodynamic fields. This 
is due to collisional cooling, which sets a minimum value for the relative size of the steady 
state gradients. An example of this situation corresponds to the well-known (steady) 
simple shear flow state ( |Goldhirsch|2003| Santos, Garzo fc Dufty|20"04 1 . This type of flow 



can only occur when collisional cooling (which is fixed by the mechanical properties of the 
particles making up the fluid) is exactly balanced by viscous heating. Unfortunately, this 
occurs, except for the quasielastic limit, only for very extreme boundary conditions (high 
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shear). The same applies for all steady states for a granular gas heated or sheared from 
the boundaries, although the behaviour of the gradients is qualitatively different for each 



type of flow (for a case by case study, please refer to the work by Vega Reyes & Urbach 



20086 I . Another example is the case of a granular gas bounded by two opposite walls 



at rest and at the same temperature. Recent comparisons ( |Hrenya, Galvin fc Wildman 
2008 ) between molecular dynamics simulations and kinetic theory provide evidence that 
higher-order effects (i.e., beyond the NS description) can be measurable for moderate 
inelasticities. 



Nevertheless, in spite of the above cautions, the NS description is still accurate and 
appropriate for a wide class of flows. One of them corresponds to small spatial pertur- 



bations of the homogeneous cooling state for an isolated system ( Brey, Dufty, Santos & 



Kim 1998 Garzo 2005). Both molecular dynamics and Monte Carlo simulations (Brey, 



Ruiz-Montero & Cubero 1999 |Brey, Ruiz-Montero, Cubero fc Garcia-Rojo||2000 1 have 



quantitatively confirmed the dependence of the NS transport coefficients on dissipation 
(even for strong values of the coefficient of restitution) and the applicability of the NS 
hydrodynamics with these transport coefficients to describe cluster formation. Another 
interesting example consists of the application of the NS hydrodynamics from kinetic the- 
ory to characterize the symmetry breaking as well as the density and temperature profiles 
in vertical vibrated gases ( Brey, Ruiz-Montero, Moreno fc Garcfa-Rojo| 2002). In the case 
of dense gases, there is evidence of the good agreement found between the NS coefficients 
derived from the Enskog kinetic theory (iGarzo & Dufty 1999a I and computer simula- 



tions (Lutsko, Brey & Dufty 2 002| |Dahl, Hrenya, Garzo fc D ufty 2002 Lois, Lemaitre 
fc Carlson||2007| ). Similar comparisons between NS hydrodynamics and real experiments 
of supersonic flow past a wedge (Rericha, Bizon, Shattuck & Swinney 20021 and Nu- 



clear Magnetic Resonance (NMR) experiments of a system of mustard seeds vibrated 



vertically ( 


Yang, Huan, Candela, Mair & Walsworth 2002 


Huan, Yang, Candela, Mair & 


Walsworth 2004 


1 have shown both qualitative and quantitative agreement. Furthermore, 



the Navier-Stokes hydrodynamics accounts for a variety of properties of granular flows 
that are preserved beyond NS order. For instance, the curvature of the temperature pro- 
file: the same types of steady profiles predicted by NS hydrodynamics have been observed 
in DSMC simulations beyond the NS domain, as it is reported in a recent theoretical 



work where new steady states (exclusive to granular gases) have been found (Vega Reyes 
& Urbach 2008&). Therefore, the NS equations with the transport coefficients derived 



here can be considered still as an useful theory for a wide class of rapid granular flows, 
although more limited than for ordinary gases. 

The plan of the paper is as follows. In §|2] the full transport coefficients of the mixture 
are given in terms of the solutions of a set of coupled linear integral equations previously 
derived by Garzo & Dufty (20021. These integral equations are approximately solved 
by means of the modified first Sonine approximation in §[3] where explicit forms for 
the transport coefficients are provided. Technical details of the calculations carried out 
here are given in two Appendices and in a Supplementary Material submitted in the 
online version of the paper. Next, in §[3] the results obtained for these seven transport 
coefficients from the standard and modified Sonine approximations are compared for 
several cases. In addition, both theoretical approaches are also compared with available 
and new simulation data obtained from numerical solutions of the Boltzmann equation by 
using the direct simulation Monte Carlo (DSMC) method ( Bird|1994 l in the cases of the 
diffusion and shear viscosity coefficients, both for two- and three-dimensional systems. 
The paper is closed in §[5] with a discussion of the results presented in this paper. 
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2. Navier-Stokes transport coefficients for a granular binary mixture 

We consider a binary mixture composed by smooth inelastic disks (d = 2) or spheres 
(d = 3) of masses mi and m 2 , and diameters o\ and a-i- The inelasticity of collisions 
among all pairs is characterized by three independent constant coefficients of restitution 
c^n, «22i an d a>i2 = a2i, where ^ 1 is the coefficient of restitution for collisions 
between particles of species i and j. In the low-density regime, the one-particle velocity 



distribution functions fi(r,v,t) obey the set of (inelastic) Boltzmann equations (Gold- 
shtein fc Shapiro||l995| Brey, Dufty & Santos 199 7). When the hydrodyna mic gradients 
present in the system are weak, the CE method ( |Chapman fc Cowling|1970 1 conveniently 
adapted to dissipative dynamics provides a solution to the Boltzmann equation based on 
an expansion 

(2-1) 



c(l) 



where is the local version of the homogeneous cooling state (HCS) ( Garzo & Dufty 
1999 & I. Although the exact form of the distribution fjy is not known (even in the one- 



component case), an indirect information of the behavior of is given through its 

from its Maxwellian form can be 



velocity moments. In particular, the deviation of /, 
characterized by the fourth cumulant 

Ci = 2 



(0) 



f 



ni Tfd{d + 2) 
where V = v — U is the peculiar velocity, 



dv V 4 f^-l 



dot 



(0) 



is the number density of species i, 

U = 



1^ 



i=l 



dv TOiV / 



(0) 



is the mean flow velocity and p — mini + miTi2 is the total mass dens ity. 



The first-order distribution fj- 1 ' has the form (Garzo fc Dufty|||2002) 



(i) 



Ai • Vxi +B t -Vp + Ci ■ VT + V lM V k Ui 



where Xi = rii/n is the mole fraction of species i 

1 



2 

dn ^ 

i=l 



dv rriiV f 



(0) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



is the temperature and p — nT is the pressure. Here, n = rij is the total number 
density. The coefficients A^, Bi, Ci, and T>i are functions of the peculiar velocity V = 
v — U and the hydrodynamic fields. 

The knowledge of the distributions f^ allows us to determine the NS transport co- 
efficients. The analysis to obtain them has been previously worked out by Garzo and 



coworkers (Garzo fc Dufty 2002 Garzo et al. 2006 Garzo fc Montanero 2007). For the 



sake of completeness, the final results will be displayed in this Section. The mass flux 



ji , the pressure tensor P£ \ and the heat flux are given, respectively, by (Garzo fc 



Duftyl[2002 l 



AD - 
Ji — 



kl ' 

mim 2 n p p , 

DVxi D p Vp- —D VT, 

P P T 



AD _ 

J2 " 



.AD 
Ji i 



(2.7) 
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-T 2 D"V Xl - LVp - AVT. 



(2.8) 



(2.9) 



It must be noted that the contribution to first order in the gradients of the collisional 
cooling vanishes (C*- 1 -* = 0) by symmetry and there is no need to consider here this 
term ( Brey et «L||1998| . This is due to the fact that the cooling rate is a scalar and so, 
corrections to first ord er in gradients can arise only from the divergence of the velocity 
field. However, as Eq. (2.5 1 shows, there is no contribution to fr 1 ' proportional to V ■ U 
and consequently, £W — 0. This is special to the low density Boltzmann equation and 
such terms occur at higher densities ( Garzo fc Dufty|[l999a Garzo et «L||2007& I. 

The transport coefficients appearing in ( 2.7 H 2.9 1 are the diffusion coefficient D, the 
pressure diffusion coefficient D p , the thermal diffusion coefficient D' , the shear viscosity 77, 
the Dufour coefficient D", the pressure energy coefficient L, and the thermal conductivity 
A. These coefficients are defined as 



D 



D 1 



P 

dm2n 

mip 
dp 

m,\T 
dp 



dvV-Ai, 

dv V • B u 
dvV-Cx, 



1 



V = 



(d-l)(d + 2) 



E 



D" = - 



L = — 



1 2 
dT 2 E 2 

i— 1 

1 2 



2 



i=l 
1 2 

d^ 2 

i—1 



dv VV : V, 



dv V 2 V-Ai, 



dv V 2 V ■ B, 



dvV 2 V -Ci. 



(2.10) 

(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 
(2.16) 



As in the case of elastic collisions ( |Chapman fc Cowling| 1970| , the unknowns are the 
solutions of a set of coupled linear integral equations. Their explicit forms are given in 



the Supplementary Material file of this article and by Eqs. (46)-(49) of Garzo & Dufty 
d2002l. 



3. Modified First Sonine approximation 

The results presented in the above Section are still exact. However, to get the explicit 
dependence of the NS transport coefficients on the parameters of the mixture (masses, 
sizes, composition, coefficients of restitution), one has to solve the integral equations 
obeying A^, Bi, Ci, and as well as one has to know the explicit forms of /f " 1 an d of 
the HCS cooling rate (see Appendix A and the work by |Garzo fc Dufty (1999&I for 
more details) . With respect to the distribution ff^ , except in the high velocity region, 
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fP(V) = f i , M (V) 



4 



d{d + 2) 



where c, is denned by (2.2 1, V* = V/v , 9i = m i r ) i 1 J2 



jfn^ 1 , and 



!'() 



/2T(mi + m 2 ) 

77117712 



(3.1) 



(3.2) 



is a thermal speed. Further, 7, = Tj/T is the temperature ratio and /i,M is the Maxwellian 
distribution 



/i,Af(^) = n 



= Ui f - 



\2nTi 

The partial temperatures T; are defined as 

1 



d/2 



exp 



dri; 



dv m. l V 2 f 



The cumulant Ci measures the departure of from fi.M- The temperature ratios 7^ 
along with the coefficients c, and the cooling rate have been determined for inelastic 
hard spheres (d = 3) (Garzo & Dufty 1999&). These calculations have been extended 



(0) 



(3.3) 



(3.4) 



here in § Appendix |A| to an arbitrary number of dimensions. 

With respect to the functions {A*i,Bi,Ci,T>i}, as in the elastic case, the simplest ap- 
proximation consists of expanding them in a series of Sonine polynomials and consider 
only the leading terms in this expansion. Usually, for simplicity, the polynomials are 
defined with respect to a Gaussian weight factor. This alternative has been previously 



used (Garzo & Dufty 2002 Garzo & Montanero 20071 to get explicit forms for the NS 



transport coefficients. The accuracy of these theoretical predictions (based on the stan- 
dard first Sonine approximation) over a wide range of inelasticities has been confirmed 
by Monte Carlo simulations of the Boltzmann equation in the cases of the tracer diffusion 



coefficient (Garzo & Montanero 2007 20041 and the shear viscosity coefficient (Garzo 



& Montanero 2007 Montanero & Garzo 20031. However, as said in the Introduction, 



recent comparisons with computer simulations for the transport coefficients associated 
with the heat flux for a monocomponent gas have shown significant discrepancies for 



strong inelasticities (Brey & Ruiz-Montero 2004 Brey et at. 2005a Montanero et at. 



20071. This fact has motivated the search for new approaches, such as a modified first 
Sonine approximation ( Garzo et aZ.||2007c 



As already pointed out by Garzo et al. (2007c I, one of the possible sources of dis- 



crepancy between the standard first Sonine approximation and computer simulations for 
the heat flux transport coefficients could be the emergence of the existing non-Gaussian 
features of the zeroth-order distribution function ff \ Although the Maxwellian distri- 
bution fi t M is a good approximation to f[ in the region of thermal velocities relevant to 
low degree velocity moments (hydrodynamic quantities, mass flux), quantitative discrep- 
ancies between /^m and f\ are expected to be important when one evaluates higher 
degree velocity moments, such as the pressure tensor and the heat flux, specially for 
strong dissipation. However, the behavior of the first-order distribution ff~' in the stan- 



dard Sonine approximation ( Garzo fc Dufty|2002 ) is mainly governed by the Maxwellian 
distribution fi.m and not by the HCS distribution f-° '. A possible way of mitigating the 
discrepancies between the standard first Sonine approximation and simulations would 



be to incorporate more terms in the Sonine polynomial expansion ( |Garz6 fc Montanero 
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|2004[ ), but the technical difficulties to evaluate these new contributions for general binary 
mixtures discard th is method. Here, we follow the same route as in our previous work 



(iGarzo et aZ.||2007c| for monocomponent gases and take the distribution ff 1 ' instead of 



the simple Maxwellian form /^a/ as the convenient weight function. In this case, some 
care must be taken in the structure of the velocity polynomials chosen to preserve the 



20031. These conditions are given by 



solubility conditions of the CE method (Chapman & Cowling 1970 Garzo & Santos 



dvfP(v) = 0, 



(3.5) 



2 , 

mi l dvvf\ x \v) = , 



i=l 
2 



(3.6) 



(3.7) 



While the conditions (3.51 and fl3.6| yield the same velocity polynomials as in the stan- 



dard method, in order to preserve the solubility condition (3.7 1 one has to replace the 
polynomial 



s,<v< = [ x - mi v 2 



d+2 r 



Ti V 



(3- 



appearing in the standard Sonine polynomial expansion ( Garzo & Dufty 2002 1 by the 
modified polynomial 



5,-(Vi = ( \m l V' 1 



d + 2 



1 



1 



Ti V 



Si(V) 



d + 2 



(3.9) 



As will be shown later, this replacement gives rise to new contributions to the transport 
coefficients associated with the heat flux. 

The determination of the NS transport coefficients in the modified first Sonine approx- 
imation follows similar mathematical steps as the ones previously used in the standard 
first Sonine approximation. More specific technical details on these calculations are avail- 
able as a supplement to the online version of this paper or on request from the authors. 
Here, we only display the final expressions for the NS transport coefficients in terms of 
the parameter space of the problem. 



3.1. Mass flux transport coefficients 

The mass flux contains three transport coefficients: D, D pi and D'. Dimensionless forms 
are defined by 



D 



PT 



-D* 



D„ 



P 



D* 



D' 



P 



D 



(3.10) 



mim 2 ^ pvn pvo 

where vq = na'^ 1 VQ is an effective collision frequency and o\2 = (o"i +<T2)/2. The explicit 
forms arc then 
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x 2 + nx\ 
D'* = - 



C 



-V 



VD ~2 C + 2u D J 



k D >- 



9 

(3.12) 
(3.13) 



In these equations, £* = C^V^Oj A* = m i/' m 2 is the mass ratio and the expression of vd 
is given by Eq. flBlj ) of the Appendix |Bflf] 



3.2. Shear viscosity coefficient 
The shear viscosity coefficient 77 is given by 

r/ = — {x\Iy7)l + x 2 T 2 r] 2 ) , 
"0 

where the partial contributions rj* to the shear viscosity are 

27 2 (2r 2 2 - C*) - 471T12 



7i72[C* 2 ~ 2C*(m + T22) + 4(r n T 22 - t 12 t 21 ) 
27i(2rn - C*) - 47 2 r 2 i 



??2 I1I2K* 2 ~ 2C*(ni + r 22 ) + 4(tiit 22 - ti 2 t 2 i)] : 
where the collision frequencies are given by Eqs. (B3l-(B6l. 



(3.14) 

(3.15a) 
(3.156) 



3.3. Heat flux transport coefficients 
The transport coefficients D" , L, and A associated with the heat flux can be written as 
d+ 2 n 



D" = - 



2 (mi + m 2 )vo 

d+2 T 
2 (mi + m 2 )vo 

d + 2 nT 
2 (mi + m 2 )u 



xrjf 
M12 



d* 



M21 



^l7l f + ^272 £ * 



M12 



M21 



^7l AI+a ^7| A , 



(3.16) 

(3.17) 
(3-18) 

A*12 A*21 \/ i 12 M21 / 

where /jjj = 777^/(777.; + m^) and the coefficients D* , D*, and D'* are given by Eqs. 
(3.11 1— ( 3.13 1, respectively. The expressions of the (dimensionless) coefficients d*, £*, and 
A* are 

d\ = \{2 [2u 12 Y 2 - Y x {2v 22 - 3(*)] [2/12^21 - ^11^22 + 2{vn + «^a)C* - 4 C* 2 ] 



71 72 
M12 M21 

7i _ _72_ 
M12 M21 

71 _ _72_ 
Ml2 M21 



A 

+2 



dC 

dx\ 



p.T 



(Y 3 + Y 5 ) [2v 12 v 21 + 2z,f 2 - C*(7^ 22 - 6C*)] 



-2V12 



dx 1 



(r 4 + ^ 6 ) (2i/u + 2^ 22 - 7C*) 



(3.19) 



^1 = T {^ 2y 3 [2(^12^21 - ^11^22)^22 + C*(7^11^22 - 5^12^21 + 2^ 22 

-6vnC ~ 7v2iC + 6C* 2 )] + 2F 4 ^12 [2^12^21 - 2v n v 2 2 + 2C>n + ^22) - C* 2 ] 
+ 2Y 5 C [2v l2 v 21 + v 22 (2v 22 - 7C) + 6C* 2 ] - 2^**6 [2(v lx + v 22 ) - 7C*]} , (3.20) 

f The definitions of collisional frequencies and the calculation of their integrals are given in 
the accompanying Supplementary Material file. 
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A* = {-^3C* [2^12^21 + ^22(2^22 
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'22 ( 5 C* " 
2^22) 



-Y s [Av 12 v 21 (v 22 - C) + 2v - 
+V12YQ [4^12^21 + 2^ii (5C* 
Here, the Y's are defined by p7|-pl2l 
A = [4(^12^21 - ^11^22) + 6£*(^n 



7C) + 6C* 2 ] + wCY A [2(i/ii + ^22) - 7C*] 

6C* 2 )+5C* 2 (6C* 



2vn) + 2vn(7v 22 C 
hC(10^22 - 23C)]}. 



7^22)] 
(3.21) 



^22) - 9C* 2 ] [^12^21 - ^11^22 + 2C*(m 



^ 22 )-4C* J 

(3.22) 



and the (reduced) collision frequencies i/y are given by Eqs. ( B 13 1 and (B14l. The 
corresponding expressions for d 2 , ^2, A2, v 22 , and v 2 \ can be deduced from (3.19 1- |3.21 1 
by interchanging 1 <-> 2. With these results, the coefficients D" , L, and A can be explicitly 
obtained from ( 3.16 1— (3.18 1 . In ( 3.19 H3.21 1 it is understood that the coefficients D* , 
D*, and D * are given by (3.11)-(3.13), respectively. Of course, our results show that 

D is antisymmetric with respect to the change 1 <-> 2 while L and A are symmetric. 
Consequently, in the case of mechanically equivalent particles (mi = m 2 , o\ = a 2 , 
Oiij = a), the coefficient D =0. 

The expressions for the NS transport coefficients derived here for a granular binary 
mixture reduce to those previously obtained (Garzo & Montanero 2007) when one takes 
Maxwellians distributions for the reference homogeneous cooling state (c\ = c 2 = 0). 
In addition, for mechanically equivalent particles, the results obtained by Garzo et al. 
(2007c) for a single gas by using the modified Sonine method are also recovered. This 
confirms the self-consistency of the results reported in this paper. 



4. Comparison with the standard first Sonine approximation and 
with Monte Carlo simulations 

The expressions for the transport coefficients derived in the previous Section depend on 
many parameters: {xi,mi/m 2 , ai/a 2 , an, a 22l oti 2 }. Obviously, this complexity exists in 
the elastic limit as well, so that the primary new feature is the dependence of the transport 
coefficients on dissipation. Thus, to show more clearly the influence of inelasticity in 
collisions on transport, we normalize the transport coefficients with respect to their values 
in the clastic limit. Also, for simplicity, we take the simplest case of common coefficient of 
restitution (an = a 22 = ct\ 2 = a). This reduces the parameter space to four quantities: 
{xi,mi/m 2 ,<Ji/<7 2 ,a}. 

Let us start with the coefficients D, D p and D' associated with the mass flux. In Fig. 



[l]thc reduced coefficients D(a)/D(l), D p (a)/D p (l) and D *, defined by (|3_13|, are plot- 
ted as functions of the coefficient of restitution a for an equimolar mixture X\ = i of 
hard spheres (d = 3) with the same size ratio u> = a\ja 2 — 1 and three different values 
of the mass ratio /x = nii/m 2 . Here, D(l) and D p (l) refer to the elastic values of D 
and D p , respectively. The coefficient D has not been reduced with respect to its elastic 
value since D' = for a = 1. We observe that both Sonine approximations lead to quite 
identical results, except for quite extreme values of dissipation where both approaches 
present some discrepancies. This is especially important in the case /1 = 2 where the dis- 
agreement is, for instance, about 15% and 32% for the coefficient D p at a = 0.1 and 0.2, 
respectively. In order to test the accuracy of both Sonine approximations in the case of 
the diffusion coefficients, the dependence of the (reduced) diffusion coefficient D(a) / D(\) 
on a is plotted in Fig. [2] in the tracer limit (x\ — > 0) for two different cases when the 
tagged particle has the same mass density as the particles of the gas (i.e., fi = uj d ). 
The theoretical predictions of both Sonine approximations are compared with available 
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Figure 1. Plot of the reduced coefficients D(a)/D(l) (a), D p (a)/D p (l) (b), and D*(a) (c) 
as functions of the coefficient of restitution a for hard spheres with xi = 4, <j\ = <Ji and two 
different values of the mass ratio fi = m\/ni2. The solid lines correspond to the results obtained 
from the modified first Sonine approximation while the dashed lines refer to the results obtained 
from the standard first Sonine approximation. 



(Garzo & Montanero 2004 20071 and new Monte Carlo simulations (using the DSMC 
method, Bird|1994 ). Here, the tracer diffusion coefficient has been measured in computer 



simulations from the mean square displacement of a tagged particle in the HCS (Garzo 



& Montanero 2004). It is apparent that both Sonine approximations provide a general 



good agrement with simulation data. However, the standard approximation slightly over- 
estimates the diffusion coefficient at high inelasticity, this effect being corrected by the 
modified approximation. 

The shear viscosity r\ is perhaps the most widely studied transport coefficient in gran- 
ular fluids. Here, this coefficient has been measured in simulations by means of a new 
method proposed by the authors ( Montanero et q£| |2005). This method is based on the 
simple shear flow state modified by the introduction of (i) a deterministic nonconser- 
vative force (Gaussian thermostat) that compensates for the collisional cooling and (ii) 
a stochastic process. While the Gaussian external force allows the granular mixture to 
reach a Newtonian regime where the (true) Navier-Stokes shear viscosity can be identi- 
fied, the stochastic process is introduced to reproduce the conditions appearing in the 
CE solution to Navier-Stokes order. More details on this procedure can be found in Mon- 
tanero et al. (2005). The simulation data obtained from this method along with both 
Sonine approximations are presented in Fig. [3] for disks (d — 2) and spheres (d = 3). We 
have considered again mixtures constituted by particles of the same mass density. The 
simulation data corresponding to d = 2 for a ^ 0.5 were reported by Garzo & Montanero 
(120071) while those corresponding to d = 3 and d — 2 for a ^ 0.5 have been obtained 
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Figure 2. Plot of the reduced diffusion coefficient D(a) / D(l) as a function of the coefficient 
of restitution a for hard disks (top panel) and hard spheres (bottom panel) in the tracer limit 
(xi 0) when the tagged particle has the same mass density as the particles of the gas (fi = ui ). 
Three different values of the mass ratio have been considered: fj. = 1, /J, = 4 and [i = 8. The solid 
and dashed lines represent the modified and standard first Sonine approximations, respectively. 
The symbols are DSMC results obtained from the mean square displacement of the tagged 
particle. The DSMC results correspond to fi = 1 (circles), fi — 4 (triangles), and /i = 8 (squares). 




Figure 3. Plot of the reduced shear viscosity coefficient ri(a)/n(l) as a function of the coefficient 
of restitution a for hard disks (top panel) and hard spheres (bottom panel) for an equimolar 
mixture (xi = |) constituted by particles of the same mass density (/i = u d ). Three different 
values of the mass ratio have been considered: fj, = 1, fi = 4 and [i = 8. The solid and dashed 
lines represent the modified and standard first Sonine approximations, respectively. The DSMC 
results correspond to fi, = 1 (circles), /i = 4 (triangles), and fi — 8 (squares). 
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Figure 4. Plot of the reduced thermal conductivity coefficient A(a)/A(l) for hard spheres with 
xi = j , <Ti = (72 and three different values of the mass ratio fi: fi — 1 (a), [i = 2 (b), and fx = 8 (c). 
The solid lines correspond to the results obtained from the modified first Sonine approximation, 
the dashed lines refer to the results obtained from the standard first Sonine approximation 
while the dotted lines are the results derived by neglecting non-Gaussian corrections to the HCS 
distributions. The symbols in the p anel (a) are DSMC results o btained from the Green-Kubo 
relations for a monocomponent gas (Brey &: Ruiz-Montero|2004 1. 



in this work. As in the case of a single gas (Garzo et al. 2007c 
a 



we observe that up to 



0.6 the simulation data agree quite well with both theories. On the other hand, 
for higher inelasticities and in the physical three-dimensional case, the standard first 
Sonine approximation overestimates the shear viscosity while the modified first Sonine 
approximation compares well with computer simulations, even for low values of a. We 
also observe that the improvement of the modified Sonine method over the standard one 
is less clear for hard disks (d = 2) since the numerical data lie systematically between 
both theoretical approaches. However, even in this case the simulation results are closer 
to the modified ones than the standard ones. Therefore, according to the comparison 
carried out at the level of the coefficients D and rj, we can conclude that while the stan- 
dard Sonine approximation does quite good a job for not strong values of dissipation, it 
is fair to say that the modified Sonine approximation is still better (especially for hard 
spheres) since is able to agree well with computer simulations in the full range of values 
of dissipation explored. 

Let us consider finally the heat flux. As said in the Introduction, recent studies for a 



monocomponent gas (Brey & Ruiz-Montero 2004 Brey et al. 2005a Montanero et al. 



20071 have shown that the standard first Sonine approximation dramatically overesti- 



mates the a-dependence of the transport coefficients associated with the heat flux for 
high dissipation (a < 0.7). This is the main reason why new alternative approximations 




Figure 5. Plot of the reduced pressure energy coefficient L(a)n/ A(l) for hard spheres with 
xi = h, ffi = <72 and three different values of the mass ratio fi: fi — 1 (a), fi = 2 (b), and fx = 8 (c). 
The solid lines correspond to the results obtained from the modified first Sonine approximation, 
the dashed lines refer to the results obtained from the standard first Sonine approximation 
while the dotted lines are the results derived by neglecting non-Gaussian corrections to the HCS 
distributions. The symbols in the p anel (a) are DSMC results o btained from the Green-Kubo 
relations for a monocomponent gas ( Brey &: Ruiz-Montero[ 2004 ) , 



( Noskowicz et aZ.|[2007 1 have been proposed. However, in contrast to the single gas case, 
the lack of available simulation data for granular mixtures for the coefficients D" , L and 
A precludes a comparison between both Sonine approximations and computer simula- 
tions. Figures |4j [5] and [6] show the a-dependence of the reduced transport coefficients 
A(a)/A(l), L(a)n/\(1), and D" (a) / D" (1) , respectively, for d = 3, x\ — \, u = 1 and 
different values of the mass ratio. Also for comparison we show the theoretical results 
obtained by neglecting non-Gaussian corrections to the HCS distributions ff^ [i.e., the 
cumulants given in (2.2 1 arc neglected: Cj = 0]. Simulation results reported by Brey & 
Ruiz-Montero (2004) for a monocomponent granular gas are also included in the case 
H = 1, showing that the modified Sonine appoximation agrees with simulation data sig- 
nificantly better than the standard Sonine approximation, especially in the case of the 
pressure energy coefficient L. As expected, the standard and modified Sonine theories 
differ significantly for strong inelasticity, especially for mechanically different particles 
(/i 7^ 1). These discrepancies clearly justify the use of the modified Sonine approximation 
instead of the standard one to compute the dependence of the heat flux transport coeffi- 
cients on dissipation. In addition, the standard approach leads to unphysical (negative) 
values for the thermal conductivity coefficient A for low values of a, which is corrected 
by our modified Sonine approximation. This is of relevance, since it helps to clarify that 
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Figure 6. Plot of the reduced Duffour coefficient D" (a)/ D(l) for hard spheres with x± = |, 
o\ — (r 2 and three different values of the mass ratio fi: fj, = 2 (a), fj, = 4 (b), and fi = 8 (c). 
The solid lines correspond to the results obtained from the modified first Sonine approximation, 
the dashed lines refer to the results obtained from the standard first Sonine approximation 
while the dotted lines are the results derived by neglecting non-Gaussian corrections to the HCS 
distributions. 



the origin of unphysical values for the thermal conductivity at very high inelasticity is 
due to the mathematical accuracy of the Sonine approximation used and not to a more 
fundamental reason; an inherent lack of scale separation, for instance. Note also that 
the simple theoretical results ( Garzo fc Montanero||2007 1 obtained by using Maxwellian 
forms for f^ are even closer to the modified ones than those predicted by the standard 
approximation. The good agreement found in general between simulation (which yield 
the " actual" values of the transport coefficients) and the modified Sonine approximation 
is a first evidence that the Chapman-Enskog method may be safely extended down to 
very low values of the coefficient of normal restitution a. 



5. Discussion 



A modified Sonine approximation recently proposed flGarzo et al. 2007c ) for mono- 



component systems has been extended to granular mixtures in this paper. This theory 



has been mainly motivated by the disagreement found at high dissipation (Brey & Ruiz- 



Montero 2004 Brey et al. 2005a Montanero et al. 20071 between the simulation data 



for the heat flux transport coefficients and the expressions derived from the standard 
first Sonine approximation for a single gas. As we have shown, important discrepancies 
between computer simulations and the standard Sonine approximation appear only in 
the region of strong inelasticity (a < 0.7). Thus, it could be argued that the search 
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for new theoretical methods may admittedly be more a mathematical than a physical 
endeavor. However, it is still physically relevant to propose methods (Noskowicz et al. 
2007 Garzo et aL||2007c ) that produce accurate results for the transport coefficients in 
the complete range of possible values of the coefficients of restitution. As in the case of 
a monocomponent gas ( Garzo et «?.||20"07c[ ), in the modified Sonine approximation the 
weight function for the unknowns Ad, 3i, Ci 7 and T> t defining the first order distribution 
is not longer the Maxwell-Boltzmann distribution f^M but the HCS distribution 



fj°\ Moreover, in order to preserve the solubility conditions (|3.5|-(|3.7|, the polynomial 



Si(V), defined by (3.8 1, appearing in the standard Sonine polyn omial expansion must 
be replaced by the modified polynomial Si(V), defined by (3.9 1. The idea behind the 
modified method is that the deviation of /■ from /^a/ has an important influence on 
the NS distribution £ . In this context, it is expected that the rate of convergence of 

the Sonine polynomial expansion is accelerated when fj ' rather than ji_u is used as 
weight function. 

The extension of the modified Sonine approach to mixtures is not easy at all since it 
involves the computation of new complex collision integrals. However, the structure of 
the NS transport coefficients is quite similar in the standard and modified approxima- 
tions (IGarzo & Dufty||2002|l since the distinction between both approximations occurs 



essentially in the a-dependence of the characteristic collision frequencies of the transport 
coefficients vd, ry, and i/y. The forms of the NS coefficients are given by equations 
( 3. 10 1— ( 3.13 1 for the mass flux transport coefficients D, D p , and D', equations (3.14|- 



( 3. 15 6 1 for the shear viscosity rj and equations ( 3. 16 1— ( 3.21 1 for the heat flux transport 
coefficients A, L, and D". All the above expressions have the power to be explicit, namely, 
they are explicitly given in terms of the parameters of the mixture. A Mathematica code 
providing the NS transport coefficients under arbitrary values of composition, masses, 
sizes, and coefficients of restitution can be downloaded from our website [f] The fact that 
our theory does not involve numerical solutions allows us to evaluate the transport coef- 
ficients within very short computing times. For instance, the code using our theoretical 
expressions for the NS transport coefficients takes just of the order of 5 seconds in a 
standard personal computer to produce a graph similar to those in Fig. [6] of the most 
complicated heat flux transport coefficients [f] These results contrast with the method 
devised by Noskowicz et al. (2007) for a mocomponcnt gas where a system of algebraic 
equations must be numerically solved by employing the power of symbolic processors. 
This is perhaps another new added value of the method developed here for granular 
mixtures. 

In order to check the accuracy of the modified Sonine approximation, computer sim- 
ulations based on the DSMC method have been carried out in the cases of the tracer 
diffusion coefficient D and the shear viscosity coefficient rj. Although some simulation 



data for these coefficients were previously reported by the authors ( Garzo & Montanero 
2004 Montanero et a?.|2005 l, in this paper we extend those simulations to very low values 
of the coefficient of restitution (values of a typically larger than 0.1). This allows one to 
make a careful comparison between the modified and standard approximations with com- 
puter simulations. As expected, the discrepancies between simulation and the standard 
first Sonine estimates for D and rj are partially corrected by the modified approxima- 
tion, showing again the reliability of such approach for very strong values of dissipation. 



f http:/ /www.unex.es/fisteor/vicente/granular_files.html/ 
J Calculations were performed in a PC equipped with an x86 64-bit Intel® processor and a 
RAM memory of 2 GB 
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Unfortunatelly, the lack of available simulation data for the coefficients A, L and D" 
corresponding to the heat flux precludes a comparison between theory and simulation 
for these transport coefficients, except in the single gas case ( Brey fe Ruiz-Montero|20(M 



Brey et aL|2005a" l where Figs. [5] and [Fj] show again the superiority of the modified Sonine 
approximation. We expect the present results for binary granular mixtures stimulates the 
performance of such simulations to assess the degree of accuracy of the modified Sonine 
method for the heat flux transport coefficients. 

Hydrodynamic theories based on kinetic theory tools have been successful at predicting 
not only rapid granular flows (where the role of the interstitial fluid is assumed negligible) , 
but have also been incorporated into models of high- velocity, gas-solid systems. These 
kinetic theory calculations are now standard features of commercial and research codes, 
such as Fluent® and MFIX (http://www.n1fix.0rg/l. These codes rely upon accurate 
expressions for the transport coefficients, and a first-order objective is to assure this 
accuracy from a careful treatment. As is shown in this paper, the price of this approach, 
in contrast to more phenomenological approaches, is an increasing complexity of the 
expressions for the transport coefficients as the system becomes more complex. In this 
context, we expect the theory reported here for granular binary mixtures to be quite 
useful for the above numerical codes. 

The present results can be also applied to determine the dispersion relations for the 
hydrodynamic equations linearized about the homogeneous cooling state. Some previous 



results ( Garzo et al. 2006 1 based on the standard Sonine method have shown that the 
resulting equations exhibit a long wavelength instability for three of the modes. The 
objective now is to revisit the above problem by using the a-dependence of the transport 
coefficients obtained in this paper (if this should be done, then for consistency reasons, 
the linearized Burnett corrections to the cooling rate would need to be considered in 



the hydrodynamic equations, since V • j '\ , ^iPj, an d V • are of second order 
in the gradients, following the analysis by Brey et al. (1998 1). Another possible open 
problem is to extend the present results (which are restricted to a low-density granular 
mixture) to finite densities in the framework of the Enskog kinetic theory. Given that 
the NS transport coefficients have been recently obtained from the Enskog equation 
(Garzo et al. 2007a|& I by means of the standard Sonine approximation, it would be 
interesting to compare again the above results with those derived from the modified 
first Sonine method. Finally, it would also be interesting to devise a method to measure 
by means of Monte Carlo simulations some of the NS transport coefficients of the heat 
flux in a granular binary mixture. This method could be based on the method proposed 



by Montanero et al. (2007) for a single gas where the application of a homogeneous, 



anisotropic external force produces heat flux in the absence of gradients. We plan to 
carry out such extensions in the near future. 
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Appendix A. Homogeneous cooling state 

In this Appendix, the expressions of the cooling rate Q ' an( A the fourth-degree cu- 
mulants are given. These expressions were reported by |Garz6 fe Dufty ( |1999&[ ) for 
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inelastic hard spheres (d = 3). Here, we extend these expressions to an arbitrary number 
of dimensions d. 

By using the leading Sonine approximation (3.1l for f[ ^ and neglecting nonlinear 

terms in c\ and c 2 , the (reduced) cooling rate CI = Ci°^ l v o (where v§ — naf^Vo) can be 
written as 

Cl =ClO + CllCl+Cl2C 2 , (Al) 



where 



Cio= dm Xl wJ 1 (1_aii)+ ^r7ir^ 21 



i 



1/2 



x — — (l + a 



13)02 



"1/2 



1 - 2^21(1 + a 12 )(l + i 



(A 2) 



o_.(d-l)/2 / _ \ 



16v^dT(f) W 
^(d-D/a (l+^T^ /2 



8dr (|) ^ 0V2 
[2(3 + 40)-3M 2 i(l + a 12 )(l + 0)], (A3) 



Ci2 = - ^7gy Wi (l + a 12 )e 2 1/2 [2 + 3 f i 21 (l + a 12 )(l + e)}. (A 4) 

In the above equations, puj — m, l l(m i + mj), Q\ = l/(/i2i7i), #2 = l/(Mi272), — 61/6 2, 
and v = y // 2T(nii + m 2 )/m 1 m 2 . The expression for can be easily inferred from ( A 2 1- 
(A4| by interchanging 1 and 2 and setting — > 

The coefficients Ci and c 2 are determined from the Boltzmann equations by multiplying 
them by V 4 , and integrating over the velocity. After some algebra, when only linear terms 
in ci and c 2 are retained, the result is 

_ AG -ED _ BE — AF 

Cl ~ BG-DF' ° 2 ~ BG-DF' (A5) 

where 

d(d+2) d(d+2) 
2lP — ^ 10 ~ Al ' h %p — C20~A 2 , (Ab) 

D _ d(d+2) d(d + 2) 

u - — — Ci2+Ai 2 , t — — ^ — C21+A21, (A 1) 

In the above equations, Ai, An and Ai 2 are given by 

/a 1 \ d - 1 3 + 2d+2a 2 1 

-0r 5/2 ^ 2 (i + e)- 1/2 M21 (i + a 12 ) 



r(I) 

x j-2[d + 3 + (d+2)0]+/i 21 (l + ai 2 )(l + 0) (ll + d+ - ^^ 3 +6 < 
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-24/4 (1 + a 12 f (1 + 8f (5 + 49) + 30/4 (1 + «i2) 3 (1 + ^) 3 } , 



(A 10) 



\ 1 • = !L 4^ ^ 5/2 ^ 2 e 2 (1 + ey 5/2 M21 (1 + ai 2 ) 



ier(f) 

x {2 [d - 1 + (d + 2)9} + 3/iai (1 + «ia) (1 + 0) [d - 1 + (d + 2)0] 
-24/4 (1 + a 12 ) 2 (1 + 0) 2 + 30/4 (1 + «i2) 3 (1 + 6 1 ) 3 } • 



(AH) 



As before, the expressions for A 2 , A 22 and A 21 are easily obtained from (A 9 1— (A 11 1 
by changing 1 <-> 2. In the case of a three-dimensional system (d = 3), all the above 
expressions reduce to those previously obtained for hard spheres (Garzo & Dufty 1999&) 
The dependence of the temperature ratio 7 = T\/T 2 on the parameters of the mixture 
is determined by requiring that the partial cooling rates £j for the partial temperatures 
Tj must be equal, i.e., 



(0) 



AO) 
^2 



c 



(o) 



(A 12) 



Once equation (A 12 1 is solved, the cumulants c 2 ; are explicitly obtained by substituting 
7 into (|A5l. 



Appendix B. Expressions for the collision frequencies 

The expressions of the different collision frequencies appearing in the NS transport 
coefficients are explicitly given in this Appendix. For their definitions please refer to the 
Supplementary Material file. The collision frequency vrj associated to the the mass flux 
is given by 



2vr (<i-l)/2 
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(Bl) 



It must be noted that the expression (B 1 1 for vr> has been obtained by considering only 
linear terms in c\ and c 2 . 

The expressions of the collision frequencies Ty corresponding to the shear viscosity are 
very long and so, for the sake of clarity, we have preferred to present them in terms of 
the operators A, defined as 
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Using (B2|, the collision frequencies rn and t\ 2 are given by 

d-l 



TU ~ d(d + 2)T(l) Xl {a 12 ) 
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The corresponding expressions for t 22 and r 2 i can be easily inferred from (B4|-(B6l. 

Let us consider now the NS transport coefficients d* , £* and A* of the heat flux. The 
quantities Y$ appearing in the expressions ( 3. 19 H 3.201 for these coefficients are given 
by 
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where 
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In the above equations we have introduced the quantities 

4(01,02) = (d + 2) (2/3 12 + 2 ) + M2i(0i + a ) x 

{(d + 2)(1 - aia) - [(11 + d)a 12 - 5d - 7]/3i 2 0r 1 } 



(B16) 



+3(d+ 3)/3 1 2 2 1 1 + 2/4 (2a 



4 



d + 3 



-«12 



+ ^2) 2 , 



S(0i, 2 ) = (d + 2)(2/3 12 - 00 + /x 21 (^ + 2 ) x 

{(d + 2)(l - ai 2 ) + [(11 + d)a i2 - 5d- 

d + 3 



(B 17) 



-3(d+3)/? 2 2 2 - 1 - 2^ 1 



-aia 



+ d+ 1 8 2 - 1 {d 1 + 8 2 f 



d+3 



aia + d + 1 [(d+ 2)0 X + (d + 5)0 2 ] 



E{6 X ,B 2 ) = 2^0^(01 + 2 ) 2 [2a{ 2 - 

-IHx{9i + 2 ) {/3i 2 0r 2 [(d + 2)0! + (d + 5)0 2 ][(11 + d)a 12 - 5d - 7] 
-0 2 0r 1 [2O + d(15 - 7aia) + d 2 (l - a 12 ) - 28a 12 ] - (d+ 2) 2 (1 - aia)} 
+3(d + 3)/3 1 2 2 0r 2 [(d + 2)0! + (d + 5)0 2 ] 
+2{3 12 8i 1 [(d + 2) 2 0i + (24 + lid + d 2 )0 2 ] 

+(d + 2)0 2 0^[(d + 8)0! + (d + 3)0 2 ], (B 18) 
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Here, 0x2 = A*i2#2 — At2i#i- From (B 13 |-( |B 19 1, one easily gets the expressions for w 2 i, 
^22 and V2x by interchanging 1 2. 

Finally, note that the expre ssions ([BT|, (|B3|, (|B4|, pl~3| , ( |B14| , and ( |B 15| ) reduce 
to those previously obtained (Garzo & Montanero 2007[ ) when one takes Maxwellians 
distributions for the reference homogeneous cooling state , i.e., when c\ = c 2 = 0. 
Moreover, in the case of mechanically equivalent particles, the expressions of vb, and 
Uij are consistent with those recently obtained for a monocomponent gas by using the 
modified Sonine approximation ( Garzo et (d.||2007~ ). 
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